Spin-orbit coupled Bose-Einstein condensate in a tilted optical lattice 
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Bloch oscillations appear for a particle in a weakly tilted periodic potential. The intrinsic spin 
Hall effect is an outcome of a spin-orbit coupling. We demonstrate that both these phenomena can 
be realized simultaneously in a gas of weakly interacting ultracold atoms exposed to a tilted optical 
lattice and to a set of spatially dependent light fields inducing an effective spin-orbit coupling. It is 
found that both the spin Hall as well as the Bloch oscillation effects may coexist, showing, however, 
a strong correlation between the two. These correlations are manifested as a transverse spin current 
oscillating in-phase with the Bloch oscillations. On top of the oscillations originating from the 
periodicity of the model, a trembling motion is found which is believed to be atomic Zitterbewegung. 
It is argued that damping of these Zitterbewegung oscillations may to a large extent be prevented 
in the present setup considering a periodic optical lattice potential. 

PACS numbers: 03.75.Mn, 67.85.Hj, 71.70.Ej 



I. INTRODUCTION 

Recent years have seen an increased interest in many- 
particle spin dynamics, also termed spintronics [l[. Most 
of this research have been devoted to condensed mat- 
ter systems, where the electron spin constitutes the spin 
degree of freedom [2j. However, related or equivalent 
phenomena may occur by considering the polarization of 
light [|| or internal Zeeman levels of cold atoms exposed 
to optical light fields [4j or spatially varying magnetic 
fields @. 

In many of the spintronics schemes developed in the 
condensed matter community, a spin-orbit (SO) coupling 
plays a crucial role for the dynamics. The particular 
form of the SO coupling can be of either Rashba [6| or 
Dresselhaus Q type, both being frequently analyzed in 
terms of an effective gauge force In other words, 

one can interpret the SO couplings as giving rise to an 
effective spin-dependent Lorentz-type force. Thus, the 
SO coupling can serve as a spin-filter or a Stcrn-Gerlach 
apparatus. As the manufacturing techniques improves, 
undesired effects deriving from crystal impurities can 
be greatly suppressed and several successful experiments 
demonstrating the SO coupling and its outings have been 
reported recently Ultracold atoms in optical lat- 

tices [ill) EH possess another possibility to implement 
effective SO couplings, and the study of it within a clean 
environment. While in condensed matter settings where 



'Electronic address: jolarson@fysik.su.se 



the coupling is an intrinsic property of the systems, for 
cold atoms it has to be externally created via spatially 
varying laser fields. 

Ever since the pioneering experiments on Bose- 
Einstein condensation (BEC) [ijj, the field has witnessed 
a tremendous progress both experimentally and theoret- 
ically [l4|. The dispersive dipole coupling between light 
and atoms enables the possibilities to study matter waves 
in periodic potentials and the experimental realization 
of model systems discussed within theoretical solid state 
physics [ll(. Today, both the weakly interacting atomic 
gases as well as more strongly interacting ones can be re- 
alized experimentally. Furthermore, utilizing the internal 
atomic hyperfine levels enable studies of the interplay be- 
tween internal and external degrees of freedom 15]. One 
interesting feature is the realization of artificial gauge po- 
tentials by means of coupling of internal atomic Zeeman 
levels [l6l4l9| . The technique of light induced gauge po- 
tentials allows for non-Abelian behavior to be analyzed 
and the SO couplings are one such example 0, HJ . 

In this paper, we investigate the dynamics of a spinor 
BEC subjected both to a two-dimensional weakly tilted 
optical square lattice as well to a set of spatially vary- 
ing light fields giving rise to an effective SO coupling 
between the internal atomic states. Thereby, its dynam- 
ics will be driven by two mechanisms, namely; (i) Bloch 
oscillations (BOs) arising due to the tilting of the lat- 
tice, and (ii) an intrinsic spin Hall effect (SHE) deriv- 
ing from the effective SO coupling. The two phenomena 
may be analyzed in terms of an extrinsic force along the 
longitudinal direction of the lattice tilt and an intrinsic 
transverse force. Being orthogonal to each other, these 
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forces render a complex correlated motion of the BEC in 
the lattice potential. Furthermore, the particular form 
of the SO coupling renders Dirac cones characterized by 
linear dispersions reminiscent of relativistic electrons. As 
an outcome of these, the evolution possesses as well a 
trembling Zitterbewegung motion p2l | . Normally, such ef- 
fect damps out very fast [23| , but we show that in a lat- 
tice situations as in the present paper it may indeed be 
long-lived and hence more easily observable. Moreover, 
the optical lattice can serve as trapping potential for the 
atoms in the BO regime and no confining Harmonic trap 
is needed which would typically prevent buildup of spin 
currents [24|. 

While SO couplings arise as an intrinsic property in 
many crystals, in cold atom systems it has to be imple- 
mented externally. This, at the one hand may be a source 
for experimental errors, but at the other it gives a handle 
for the SO coupling; its strength and form can be tuned 
experimentally. The same holds for the optical lattice, 
the lattice strength and geometry are both control pa- 
rameters, which has made it possible to observe BOs for 
BECs in optical lattices (25|. Due to the large atomic 
mass compared to the mass of the electron, Zitterbewe- 
gung is expectedly more clearly manifested in atomic sys- 
tems. The first experiments detecting Zitterbewegung in 
atomic systems were recently carried out [23], [2(| . A dif- 
ferent but related setup can be found in Ref. @, where 
a spinor condensate accelerating in a periodic magnetic 
field was studied. We note that in Ref. the spatially 
dependent magnetic field gives rise to an effective gauge 
force. However, Bloch oscillations were not considered, 
and moreover, in the present work the periodic potential 
is external and independent of the SO coupling. 

The outline of the paper is as follows. In Sec. HH we 
introduce the model system and outline how effective SO 
couplings are rendered in such systems. To start with 
we analyze the free case in which the problem is analyt- 
ically solvable and demonstrate the SHE and the Zitter- 
bewegung effect. Applying the acceleration theorem we 
also solve the problem under constant acceleration. The 
dispersions, once an optical lattice has been included, 
are solved numerically displaying the Dirac cones. Sec- 
tion Mil presents the main analysis. The full numerical 
results are presented, where we demonstrate the presence 
of a SHE and how it is modified by the longitudinal os- 
cillatory motion. In particular, it is shown that the Hall 
current will also oscillate, however with an increasing os- 
cillation amplitude. On top of the oscillations deriving 
from the motion in a periodic potential, we find a per- 
sisting Zitterbewegung effect. Finally, we conclude with 
a summary in Sec. IIVI 



II. MODEL SYSTEM 

Rashba @ and Dresselhaus Q SO couplings ap pea r 
in a variety of solids, ranging from semiconductors [271 ] . 
thin magnetic films 128], to graphene [29| . The SO cou- 



plings bring about a degeneracy in the energy dispersions. 
Around these degeneracies, the dispersions are linear, 
reminiscent to that of relativistic particles, and there- 
fore they are called Dirac cones. The relativistic prop- 
erties have been confirmed experimentally in graphene 
either by angle resolved photoemission spectroscopy or 
by considering the conductivity and cyclotron motion 
of electrons [3(| The SO coupling and its corre- 

sponding degeneracy induces as well a non-trivial effec- 
tive gauge structure, giving rise to phenomena such as 
the intrinsic spin and anomalous Hall effects [1, H|. In 
the field of cold atoms, the first light induced effective 
magnetic fields were demonstrated in Ref. [l{|, and sig- 
natures of non-Abelian structures were recently observed, 
see [32| and reference therein. BOs with ultracold atoms 
have already been realized experimentally by numerous 
groups (25l . l33j . We proceed by presenting a scheme which 
realizes the desired SO coupling in an atomic tripod sys- 
tem. 



A. Light induced effective spin-orbit coupling 

The idea of effective gauge potentials arising from adia- 
batic evolution dates back to the seminal works by Mead 
on vibrations of molecules (34J. Berry [35j generalized 
this finding by demonstrating the existence of an effective 
gauge structure in the context a quantum adiabatic evo- 
lution. Wilczek and Zee [36| further generalized Berry's 
result and found an effective non-Abelian gauge poten- 
tial accompanying adiabatic evolution of degenerate en- 
ergy eigenstates. Dum and Olshanii [l6[ showed that 
the effective gauge potentials can be generated in sys- 
tems of ultracold atoms subjected to spatially dependent 
light fields that couple internal atomic states. Follow- 
ing an approach for creating non-Abelian gauge poten- 
tials in time-dependent systems [37) . Ruseckas et al. [TtJ 
outlined a scheme how to achieve effective non-Abelian 
gauge potentials by considering an atomic tripod system 
in spatially dependent light fields. The coupled four-level 
system possesses two degenerate dark states that, within 
certain parameter regimes, can be decoupled from the 
remaining two bright states via standard adiabatic elim- 
ination. The flexibility of this tripod system, deriving 
from different laser configurations, allows for a variety 
of phenomena to be studied, such as, e.g., relativistic 
properties [H [H HI, SHEs [1|, novel phases of the 
BEC 4l|, and topological defects as effective magnetic 
monopoles [42| or vortices [32j . 

In the present paper we consider a tripod setup sim- 
ilar to the one described in Refs. [2(J Hy. The details 
of atomic and laser configurations rendering the gauge 
potentials can be seen in Fig. [TJ The internal dynam- 
ics, is driven by three lasers that couple the three lower 
atomic states |1), |2), and |3) to the excited state |0). 
Parameterizing the laser fields as Oi = O sin 8 cos ipe lS1 , 
02 = sin 9 sin tpe lS2 , and 03 = cos 9e lSs , such that 

o = 



Oi| + IO2I + I O3 1 , the two degenerate dark 
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(b) 



FIG. 1: Atomic (a) and laser (b) configurations. The specific 
laser arrangement, made up of standing and traveling waves, 
is described in detail in the text. 



states of the corresponding Hamiltonian can be written 
as 
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where Ski = Sk ~~ Si. Provided the field amplitude f2 
is much larger than the atom-light detuning A (as de- 
picted in Fig. [1]), the two dark states can be safely elimi- 
nated from the two bright states, resulting in an effective 
two- level system [llj ■ To derive an effective SO coupling 
between the two internal dark states and \u2), we 
choose Si = S 2 , 5*3i = VsV, <P — v v x, and 9 <E [0,7r]. 
This can be achieved by coupling the two ground states 
|1) and |2) to the excited state |0) with two lasers in- 
tersecting with an angle £ and propagating in the xy- 
plane. These lasers correspond to a standing wave in the 
x-direction and a plane wave in the y-direction. A third 
laser, couples the states |3) and |0), and is a plane wave 
in the yz-plane with an angle £ of propagation relative to 
the j/-axis. In terms of the angles £ and £, and the wave 
numbers kj = \kj\, j = 1, 2, 3, we have 



S ± = S 2 = fcycos(C/2), 
S3 = hycos^, 
tp = 2fc£sin(C/2), 



(2) 



with k = k\ — k 2 . Performing the adiabatic elimina- 
tion, taking into account that f2 ^> |A| and the fact 



that we consider ultracold atoms so that excitations due 
to atomic motion can be neglected, the effective single 
particle two-level Hamiltonian becomes in dimensionless 
variables [20? 



H 



(0). 



'cff — 

where So 
(v v /2) cos 



pI + Pi + Pi + do&z - v a x p x - v^py, (3) 



is an effective Zccman splitting, vq = 
' "''•> The Pauli matri- 



ces a 



.1 • 



J 



vi = (v s /2) sin^/2. 
x,y,z, operate in the space spanned by 
& x = |ui)(ti2| + \u 2 )(ui\, ct v = 
and a z = - \u 2 ){u2\- The 

p-dependence of the last two terms originates from non- 
adiabatic couplings 



the dark states, i.e. 
-i\v,i}(u 2 \ +i\u 2 )(ui 



431 ] . We note that these do not in- 



clude p z and as a consequence the z-dependence can be 
fully factored out. In Eq. (O, we have used dimensionless 
parameters. Explicitly, when we further on introduce a 
lattice potential whose wave number is k, we let the char- 
acteristic length scale be I = 2irk~ 1 and we let the charac- 
teristic energy scale be the recoil energy Er = h 2 k 2 /2m. 
Consequently, time t is measured in units of r = Ti/ Er. 
The superscript (0) indicates the free Hamiltonian in ab- 
sence of a spatially dependent optical lattice potentials. 

The Zeeman shift Sq splits the degeneracy of the two 
internal states and together with the SO coupling it gives 
rise to a transverse current, an intrinsic anomalous Hall 
effect. We will, however, focus on the intrinsic SHE ob- 
tained when Sq — 0. In this case, the transverse par- 
ticle current vanishes but, on the other hand, a trans- 
verse spin current persists. The So = condition can 
be achieved by internal Stark shift of the individual bare 
atomic levels [2(| |44|. Thus, in the following we will as- 
sume So — 0, and furthermore we will also consider the 
symmetric SO coupling characterized by vq = v\ = v so . 
It should be pointed out, however, that even by lifting 
these constrains an interplay between transverse and lon- 
gitudinal atomic motion will still occur. 

One important prerequisite for successful application 
of the current scheme is that the system can be prepared 
in the dark subspace spanned by {|wi), ^2)}. This issue 
has been discussed in [451 ] . The idea is to slowly turn on 
the external lasers f^, k = 1, 2, 3, such that the state fol- 
lows adiabatically states in the dark subspace, preparing 
a superposition of \u 2 )} but no bright state com- 

ponents. 



B. Analytical results without a lattice potential 

We note that we may rewrite the Hamiltonian in 
Eq. © (S = 0) as 

= (p x -A x ) 2 + (p y -A y ) 2 + $ (4) 



with the vector potential 

(Ax, Ay) = -v so (a x ,a y )/2 
and the scalar potential 

<l = --v 2 1 
2 so 



(5) 
(0) 
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These potentials constitute vector and scalar parts of 
the effective gauge potential. Although being constant 
in space, the vector potential gives rise to a non-trivial 
effective gauge field due to its non-Abelian structure, 
[Ar,.<4y] 7^ 0. Explicitly, this effective field is given by 
the antisymmetric field tensor 

Fki = dkAi - diA k - i[Ak,Ai], k,l = x, y, (7) 

resulting in P xy — \v 2 so a z and = otherwise. As a 
result, the spin-orbital motion of the electrons can be un- 
derstood in terms of an effective Lorentz-type force, akin 
to that found by Wong [46( for a color-charged classical 
particle moving in a non-Abelian gauge field. The corre- 
sponding force operator may be derived from Heisenberg 
equations of motion according to (47j 



with 



(8) 



where it is understood that the operators on the right- 
hand side are all given in the Heisenberg representation 
and the prefactor ^ in front of the acceleration is the mass 
in the chosen units. Note that the force in the y-direction 
is proportional to — p x and the force in the x-direction 
is proportional to p y . This demonstrates the transverse 
character of Fso and its equivalence to a Lorentz or Cori- 
olis force. The crucial aspect of Eq. (j8]) is, however, that 
the force is state dependent: atoms in dark state \u\) 
feel a force with opposite sign compared to those in dark 
state |u2). Therefore, the motion of the particles in the 
xy-plane may separate the two spin states implying that 
a non-zero spin current is generated. In the following, we 
shall call the x-direction longitudinal and the y-direction 
transverse. 

The effect of the above Lorentz-type gauge force can 
be examined by computing how the initial (t = 0) state 

^{Px,Py,0)=1pl(Px,Py,0)\ui) +1p2(Px,Py,0)\u 2 ) (9) 

evolves in time under influence of H^g . We assume that 
the internal and momentum states of the atom are ini- 
tially uncorrelated and that the initial momentum states 
are Gaussian wave packets of the form 



1p2(Px,P V ,0) = ^(Px,Py) 

(P^-PO) 2 +Py 

2A' 2 

e " 



1 



(10) 



where A p sets the momentum uncertainty and po is the 
average momentum along the x-direction. This choice 
corresponds to the pure internal atomic state |— ) = 
— \u2))/v2. Given this initial state, the time evolved 
state is readily obtained as 

^(Px,Py,t) = l/)l{Px,Py,t)\ul) + ^2(Px,Py,t)\u 2 Xll) 



^i{Px-,P y ,t) = —7= [cos(Ai) + isin(Ai)e l ' 
x^(p x , Py )e-^+P>, 

i>2(p x ,Py,t) = -—7= [cos(Ai) + isin(A<)e J 
Y.^{p x ,p y )e- l{ rt + v>. 



(12) 



Here, we have used the energy dispersions E± and the en- 
ergy eigenstates <fi± of the Hamiltonian in Eq. ([3]) , which 
read 

E±(Px,Py) = Pl+Pl± A, 

\<t>±) = ^(±|t) + e ie |i», (13) 



where A = v so ^jp 2 x + and tan(0) = p y /p x . 
In Fig. [2] we display the dimensionless averages 



(y)i = 2t 



d 

IP* (Px,Py,t) 7^- Ipi [Px ,Py,t)dp x dp y , (14) 



with i = 1, 2. The factor 2 derives from normaliza- 
tion of the individual wave packet components. Since 
the initial velocity is along the longitudinal x-direction, 
the deviation from (y), = is a direct consequence of 
the state-dependent Lorentz-type gauge force Fso ■ Note 
that (y)i = —(y) 2 giving a vanishing transverse particle 
current. This symmetry is broken by keeping the Zee- 
man term 8q<j z in the Hamiltonian of Eq. ([3]). The par- 
ticle current induced by such a <7 z -term goes under the 
name intrinsic anomalous Hall effect. After t ss 20, the 
separation of the two internal dark states has reached 
an asymptotic value. Until t rj 20, the transverse spin 
current shows an oscillating behavior deriving from the 
same mechanism as the one underlying relativistic Zittcr- 
bewegung in the Dirac equation [22j. In the Dirac equa- 
tion the trembling motion describes coherent coupling 
between positive and negative energy solutions, while in 
the present situation it describes population swapping 
between the internal states. In solids, the amplitude 
of these oscillations become extremely small due to the 
lightness of the electron, while in atomic or ionic systems 
the heavier masses imply larger amplitudes and possi- 
ble detection [23|, [2(|. Despite heavier masses, the rapid 
damping of the trembling motion as depicted in Fig. 1 
is a limitation also in atomic/ionic Zitterbwegung exper- 
iments [23l |2(| . In the next subsection we will argue how 
this decay may be prevented by studying the effect in the 
presence of an optical lattice. 

In order to generate a current in the system one ap- 
plies a constant longitudinal force F ex e x , e.g., induced by 
gravitation in atomic systems. In the particular case of 
solids, the force derives from a voltage across the solid, 
and the electrons reach a steady state motion in which 
the applied force is balanced by scattering against crystal 
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t 



FIG. 2: (Color online) Time evolution of {y)i,2 (green curve 
correspond to the dark state |lti) and the red curve corre- 
spond to the dark state |«2)) as defined in Eq. (|14p . The 
dimensionless parameters in Eqs. Q and (|10[) are v so = 1 
(= vi = Da), So = 0, po = 1, and A p = 0.1. The initial state 
is thus a Gaussian with a non-zero average momentum along 
the a;-direction. 



impurities. On the other hand, a system of cold atoms 
is, to a good approximation, free from impurities and one 
would therefore achieve constant acceleration. Assuming 
that the force is weak and that the total time of evolution 
is relatively short, we can perform an adiabatic approx- 
imation, consisting in replacing p x by p x — F ex t in Eq. 
(|12p . Hence, the momentum grows linearly in time in 
accordance with classical motion. How this comes about 
can be easily understood by considering the time evolu- 
tion of the Hamiltonian H = Hq — F ex x for small forces 
F ex and times t. In this regime, the time-evolution oper- 
ator can be approximated accordingly 

(j(t) = e- i6t « e iF °* M e- iAot , (15) 

so that the first exponential gives the coupled time- 
evolution while the second simply boosts the x- 
momentum by —F ex t. For longer time periods, the large 
velocity rendered by the acceleration will typically cause 
non-adiabatic excitations, which we neglect for the mo- 
ment. Similarly, for large values of F ex t, the factoriza- 
tion of the exponentials is no longer justified. The re- 
sulting averages (y)i calculated within this approxima- 
tion are depicted in Fig. [3] As before, we encounter a 
build-up of a transverse spin current but in this case the 
current changes sign after some time and finally there is 
no net spin separation. By comparison with a full nu- 
merical integration of the Hamiltonian we find that the 
imposed adiabatic approximation breaks down after rela- 
tively short times. Nevertheless, numerically we have ver- 
ified that the dynamics beyond this approximation goes 
as well asymptotically to zero for large times. We believe 
that the origin of this result is a complicated interplay 
between the transverse and longitudinal forces, together 



with the coupling between the two internal states. For 
spinlcss particles, the two forces render a spirally trans- 
verse motion and hence a net transverse current. For par- 
ticles with spin-dependent Lorentz forces F$o one would 
expect a net transverse spin current, but the coupled na- 
ture of the two spin states seems to destroy such spin 
current in the case of constant acceleration. 
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FIG. 3: (Color online) Same as Fig. [5] but with a constant 
weak force of magnitude F ex in the a;-direction added to the 
Hamiltonian ((3} and assuming the validity of the adiabatic 
approximation. The dimensionless parameters are v so = 1, 
So = 0, po = 0, A p = 0.1, and F ex = 0.01. 



C. Band spectrum and the acceleration theorem 

In the following section we numerically investigate the 
dynamics in the presence of an optical lattice. Insight 
into the evolution behavior is gained by considering the 
spectrum of the Hamiltonian including a lattice potential 
while disregarding the linear force potential. Adding to 
the SO Hamiltonian in Eq. (j3|) a spatially oscillating 
periodic potential, yields 

H = P 2 x +P 2 y + V so {PxO-x + VyOy) 

+^cos 2 (ie) + Fcos 2 (y). (16) 

The Hamiltonian is periodic and consequently possesses 
a band-gap spectrum E v (k x ,k v ) with associated Bloch 
eigenstates u^,k x ,k y (x,y), characterized by a discrete 
band-index v = 1, 2, 3, ... and quasi-momentum num- 
bers (k Xl k y ) restricted to the first Brillouin zone, i.e., 
— 1 < k x < 1 and — 1 < k y < 1. A detailed study of Bloch 
theory in the presence of SO couplings can be found in 
Ref. [43|. The first two bands E\{k x ,k y ) and E2(k x ,k y ), 
obtained from numerical diagonalization of the Hamil- 
tonian, are shown in Fig. 01 A similar band structure 
has been obtained in Ref. [39j | for cold atoms in a square 
optical lattice. We directly see the effect of the SO cou- 
pling in terms of Dirac cones at (k x , k y ) = (n, m), where 
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n, m =, — 1, 0, 1. It is further seen that the ground state 
is not the Ui^^(x, y) Bloch state, as for a regular square 
lattice lacking SO couplings. This fact is at the heart of 
the Jahn- Teller effect [49| , where the Dirac cone (or coni- 
cal intersection as it is usually referred to when it appears 
in position space rather than in momentum space) in- 
duces a symmetry breaking that lowers the ground state 
energy. Furthermore, we note that the two lowest bands 
are not gapped, while there is a gap between the second 
and third band. Indeed, due to the SO coupled two-level 
structure of the Hamiltonian, the bands come in pairs 
where each pair possesses Dirac cones for quasi-momcnta 
(k x ,k y ) — (n,m). The effective number of Dirac cones 
is four (easily seen by shifting the Brillouin zone by \ in 
both directions), compared to two Dirac cones encoun- 
tered in honeycomb lattices, like in graphene. We notice 
that recently it was shown how more complex lattice po- 
tentials may render single Dirac cones within the first 
Brillouin zone 50] . 
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FIG. 4: The two lowest Bloch bands E v (k x ,k y ) for a square 
lattice with Rashba spin-orbit coupling. In this example, 
v so — 1 and V = 5. 

In the adiabatic regime, a non-zero extrinsic force 
F ex e x is treated by using the same approximation as 
above, yielding 

u v ,k x ,k y (x,y) -^u v , kx -F ex t,k y (x,y). (17) 

In the presence of a periodic potential, this adiabatic 
assum ptio n has become known as the acceleration the- 
orem 51]. It implies that no population is trans- 
ferred between the individual bands v. As a direct 
consequence, an initial state with well localized quasi- 
momentum within one single band will show an oscilla- 
tory motion rather than a constantly accelerating one, 
once the extrinsic force F ex e x has been turned on. These 
are the well studied Bloch oscillations [52|, which have 



been experimentally verified in a variety of systems, such 
as cold atomic gases [33|, BECs I25|, semiconductor su- 
perlattices [53| , waveguide arrays |54j , and photonic crys- 
tals [11]. 

BOs are clearly an outcome of the periodic potential 
and its characteristic energy spectrum possessing gaps of 
forbidden energies at the center and edges of the Brillouin 
zone. As the amplitude of the potential V vanishes, the 
gaps close and the acceleration theorem breaks down. In 
particular, for weak lattice amplitudes, the gap between 
the first two energy bands scales as ~ V. By linearizing 
the dispersions around the band gap and replacing k x 
by k x — F ex t, one obtains a realization of the celebrated 
Landau-Zener model [56j]. Here, k x is the initial quasi- 
momentum in the ^-direction. This model is analytically 
solvable, giving the probability of population transfer be- 
tween the two states of the corresponding bands. It has 
been shown that the decay of BOs can be correctly esti- 
mated using the Landau-Zener theory [2^, [HtJ ■ 

In the present model, the situation becomes qualita- 
tively different. At the Dirac points, the gap vanishes 
and the acceleration theorem cannot be imposed. How- 
ever, from Fig. @] it is seen that the Dirac cones appear 
between the first two bands, and not between the second 
and third band. Hence, assuming a state initially pre- 
pared on the lowest band one may expect predominant 
population transfer between the v = 1 and v = 2 bands 
and not to higher bands. Thereby, further acceleration is 
hindered by the gap between the second and third band 
and spin-dependent BOs may still be encountered. In 
other words, a modified acceleration theorem involving 
pair of bands could result. In fact, as already mentioned, 
the acceleration theorem is a specific example of the more 
general adiabatic theorem, and in molecular and chemical 
physics it has long been known that adiabaticity breaks 
down in the vicinity of a conical intersection (58| . In this 
context, it should be kept in mind that adiabaticity is 
a widespread concept in physics, and the general defini- 
tion concerns systems whose dynamics is governed by an 
explicitly time-dependent Hamiltonian. Adiabaticity, as 
it appears in the acceleration theorem, is in this sense 
more reminiscent of the Born-Oppenheimer approxima- 
tion [58| . applicable to time- independent Hamiltonians 
describing light and heavy quantal degrees of freedom 
corresponding here to the spin and the spatial motion, 
respectively, of the particles. 



III. NUMERICAL RESULTS 

Our numerical analysis will be restricted to a system 
consisting of weakly interacting atoms, for which a mean- 
field description is expected to capture the phenomena 
we are interested in. In the s-wave scattering regime, we 
thereby consider dynamics rendered by a spinor Gross- 
Pitaevskii equation 0, [2l| 
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d T , , { d 2 d 2 



\ dx 2 dy 2 

gun^x, y,t) + gi 2 n 2 (x, y, t) 




+ V (cos 2 (a;) + cos 2 (y)) - F ex x + v s 



l dx + dy U 



d a_ 



(18) 







gi2ni(x, y, t) + g 22 n 2 (x, y, t) 



>*(x,y,t). 



Here, 

®[x,y,t) 



(x,y\( Ul \*(t)) 
{x,y\(u 2 \V(t)) 





i/ji(x,y,t) 




_ ip2(x,y,t) _ 



(19) 



is a two-component spinor wave function of the dark 
states, rii(x,y,t) — \ipi(x, y, t)\ 2 , F ex the external longi- 
tudinal force, and gx\, g 22 and 512 are the effective s-wave 
scattering amplitudes related to the scattering lengths 
as gi = N4iTh 2 ai/(3Eji\/27rl z m) : where l z is the width 
of the single site wave function in the z-direction. This 
definition of gi implies that the wave function has been 
normalized as J \^(x, y, t)\ 2 dxdy = 1. In the basis of 
dark states, the scattering lengths are effective ones dif- 
ferent from the regular hyperfine scattering lengths. Nev- 
ertheless, the ratio between the two scattering amplitudes 
should approximately be the same regardless of basis, 
and, in particular, for a spinor condensate such as that 
of 87 Rb one has gu ~ g 22 ~ gi 2 [59| , We will restrict the 
analysis to repulsive interactions, gi 2 , g 22 , g\ 2 > 0, and 
we will furthermore assume gu = g 22 = g\ 2 = g [59j. 

In the absence of a lattice, the SO co upling g ives rise to 
a non-Abelian Mead-Berry curvature [3J-|36(, which re- 
sults in a transverse splitting of spin components, the 
intrinsic SHE. However, in Figs. [5] and [3] the motion 
along the longitudinal x-direction was either constant or 
monotonously accelerating, while for a weakly tilted lat- 
tice, where BOs dominate the longitudinal motion, it is 
neither clear how the transverse spin current nor how the 
trembling Zitterbewegung motion will be manifested. For 
instance, will there be sufficient time for a transverse spin 
current to be established within a single Bloch period? 
In Ref . [go| , the intrinsic anomalous Hall effect was stud- 
ied in a single particle system confined by an harmonic 
potential and hence the motion was as well oscillatory as 
in the case of BOs. The transverse current appeared as a 
rotation on top of the harmonic oscillatory motion, i.e., 
in polar coordinates (r, ip) the gauge potential causes a 
non-zero current in the (^-direction. Such current would 
indeed be zero for an SO coupling possessing an Abelian 
structure rather than non-Abelian, e.g. for couplings on 
the form v so (a a p x + a a p y ) with a = x, y, z. It is clear 
that a corresponding rotational current cannot be strictly 
encountered in the present system since exact polar sym- 
metry is broken by the tilting of the lattice and by the 
lattice in itself. 

The time evolution of the system is solved by employ- 
ing the split-operator approach (6lj . The method relies 



on factorizing the time-evolution operator into a spatial 
part and a momentum part valid in the limit of infinitely 
small time-step propagation. One convenient feature of 
this method is that the wave function is in principle ob- 
tained at any instant of time in both position ty(x,y,t) 
and momentum ^{p x ,p y ,t) space. Furthermore, for a 
system bounded from below, the approach is capable of 
giving an approximate ground state of the system by sim- 
ply propagating it sufficiently long in imaginary time. 




FIG. 5: (Color online) Absolute value of the ground state 
wave function in momentum space (a) and (c), and in position 
space (b) and (d). The dimensionless parameters are V = 5, 
Fex = 0, and g = 1 for all plots, v so — 2 in the upper plots (a) 
and (b), and v so = 1 in the lower plots (c) and (d). Thus, the 
lower plots correspond to the dispersions depicted in Fig. [4] 



Our numerical considerations will focus on propagation 
of an initial (approximated) ground state for the untilted 
(F ex — 0) lattice. That is, we consider the ground state 
of a spin-orbit coupled BEC in an optical lattice, and 
then we introduce a quantum quench by suddenly switch 
on a linear force in the x-direction. Hence, we begin by 
propagating an ansatz function in imaginary time with 
zero extrinsic force F ex = 0. When convergence of the 
imaginary time evolution has been reached, we switch 
to real time propagation, including the extrinsic force by 
adding the potential term —F ex x to the Hamiltonian. As 
predicted by the acceleration theorem, for short times the 
weak forces F ex will induce a linear increase in the quasi- 
momentum k x . In contrast to BOs in a regular square 
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lattice, here the SO coupling implies that also the quasi- 
momentum k y is indirectly affected by the force F ex . At 
the same time, the build up of motion in the y-direction 
induces an intrinsic force Fso along the longitudinal x- 
direction of BOs. In addition, the linear dispersions in- 
duce a trembling Zittcrbewegung motion. It is thereby 
likely to expect a complex motion arising from the inter- 
play between the two forces related to SHE and BOs and 
the outcomes from linear dispersions. 



Figure \5\ displays the absolute values of the initial 
ground state wave function achieved via imaginary time 
evolution for F ex = 0. The momentum wave functions 
depicted in (a) and (c) demonstrate the Jahn- Teller ef- 
fect as the ground state populates non-zero momentum 
states [49]. The lower two plots (c) and (d) correspond 
to the parameters of the two lowest Bloch bands shown 
in Fig. SJ The fact that the ground state obtained from 
this numerical method is not symmetric, as seen from 
the plots, is believed to have a two-fold cause. First, it 
seems to derive from the fact that the imaginary time 
evolution is taken only up to a point where the state 
evolution has approximately converged, and, secondly, 
from the fact that even small fluctuations can result in 
an asymmetric ground state, due the four-fold degener- 
acy of the ground state of our Hamiltonian. In fact, the 
complete convergence of the ground state has not been 
reached since the position space wave functions depicted 
in (b) and (d) are localized in contradiction to the true 
ground state having a Bloch structure, and furthermore 
since the momentum wave functions have a finite width. 
This, however, should not be seen as a shortcoming since 
normally one would cool the atomic gas in the presence 
of a confining potential which would result in a local- 
ized gas. In addition, if the propagated state extends 
over the full lattice one would run into troubles when 
quantities such as the mean spatial position is calculated 
numerically. Note that the momentum wave functions of 
(a) and (c) are different from the quasi-momentum as is 
evident in the multiple Bragg components separated by 
A.p x , Ap y — ±1. In other words, even though the mo- 
mentum wave functions extend beyond the first Brillouin 
zone, it populates, to a very good approximation, only 
the lowest band. We also note that the position space 
ground state wave functions, shown in (b) and (d), pre- 
dominantly populate the individual lattice sites. Hence, 
only judging from Figs. [5] (b) and (d) it would be dif- 
ficult to predict the presence of a SO coupling and the 
Dirac cones that derive from it, which instead manifest 
themselves in momentum space in terms of a Jahn- Teller 
effect. 



In order to analyze the interplay between the trans- 
verse and longitudinal forces, F$o and F ex e Xl we calcu- 




100 200 300 400 500 600 700 800 900 1000 
t 



FIG. 6: (Color online) Average transverse position (y) (a), 
average longitudinal position (x) (b), and average intrinsic 
force -Fso in the transverse direction (c). Red lines correspond 
to the dark state |iti) while green lines to |«2). For {£) the 
two lines are identical because of symmetry reasons. The 
dimensionless parameters are the same as in Fig. [5] (a) and 
(b), but now with a non-zero extrinsic force, F ex = 0.01. 



late 



(y)i = jjr J \A(x,y,t)\ 2 y 



dxdy, 



(®)i = Jf. I \^i(x,y,t)\ 2 xdxdy, 



(20) 



9boJ^ I \^i(Px,Py,t)\ 2 p x dp x dp yi 



with i = 1, 2 corresponding to averages in terms of 
the two dark states \ui) and \u2), respectively, Nj — 
JJ \ipi(x,y,t)\ 2 dxdy, and the last quantity is the intrin- 
sic force in the y-direction. We find that (x)i = (x)2 (up 
to numerical accuracy) and therefore plot only one of the 
two components. 

Our numerical results are presented in Figs. |5] and [7] 
The difference between the two figures lies in the scatter- 
ing amplitude g. In Fig. [7J the effective scattering am- 
plitude g is already as large as 10 recoil energies Er 
and we can thereby conclude that the spin current is 
rather insensitive to the value of g. It is generally known 
that strong non-linearity can qualitatively change system 
properties. In BECs it leads to phenomena such as soli- 
ton and vortex creation. Normally the greatest im pac t 
of the non-linearity is seen around curve crossings [62j. 
Beyond a critical value of the non-linearity (in our case 
g > gcrit), in one dimension swallowtail loops are formed 
in the vicinity of the curve crossings characterizing ad- 
ditional stationary non-linear states, e.g. solitons. Up 
to this critical value the dynamics seems rather insensi- 
tive at least on a mean-field level. However, formation 
of these loops beyond the critical value implies break- 
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down of the acceleration theorem [63|. To date, it seems 
that analysis of these effects have been restricted to one 
dimensional situations and it is therefore not fully un- 
derstood how non-linearity affects the structure of the 
Dirac cones appearing in two dimensions. Nonetheless, 
we have numerically found that the general structure of 
Figs. H] and [7J drastically change for very large values of 
g. All present BO experiments, however, are performed 
in the weakly interacting regime where it is known that 
loops are absent, and we thereby restrict our analysis to 
these experimentally more relevant parameter ranges. 

We see from both Fig. [6] and [7] that the x and y mo- 
tions are indeed correlated; when the Bloch oscillation 
reaches its turning points, the spin separation in the y- 
direction is roughly maximal and as the velocity in the 
longitudinal x-direction changes sign, the spin separation 
decreases due to the reversed sign of the intrinsic force. 
When the longitudinal velocity is the greatest (equivalent 
to maxima of the intrinsic force Fso), t Qe spin current is 
reversed. More precisely, as the atoms begin to acceler- 
ate in the tilted lattice the quasi-momentum k x increases 
linearly in time and a spin current is built up in the trans- 
verse direction. However, after half a Bloch period the 
quasi-momentum reaches the Brillouin boundary and the 
x-motion is stopped and then reversed. Consequently, 
also the spin current is reversed and one encounters a de- 
creasing spin current in the y-direction. Finally, the spin 
current has dropped to zero and builds up in the oppo- 
site direction instead. Hence, the spin current follows the 
same oscillatory behavior as the BOs. 

One may argue that a large spin separation should be 
hindered by the alternating sign of the velocity in the 
x-direction. Nonetheless, from our numerical results de- 
picted in Figs. [6] and [JJ we note that the spin separa- 
tion in the y-direction is relatively large. Indeed, (y) 
exceeds several site-spacings and furthermore this sepa- 
ration amplitude increases in time. The fact that neither 
the oscillatory spin nor the oscillatory longitudinal cur- 
rent show a clean perfect harmonic pattern derives from, 
at the one hand possible non-adiabatic effects and the 
breakdown of the acceleration theorem, and at the other 

(x) 

that the transverse current induces an intrinsic force FgQ 
in the longitudinal direction which adds to the extrin- 
sic force F ex e x . The full coupled dynamics is therefore 
rather complex, but nonetheless, our figures make clear 
the overall in-phase oscillations between the two currents. 

On top of the oscillations of (x)i, there is a slight drift 
toward negative values of x. This is an outcome of non- 
adiabaticity and breakdown of the acceleration theorem. 
However, this drift is rather small and the dynamics is 
therefore predominantly adiabatic. This drift can be sup- 
pressed by increasing the lattice amplitude V on the cost 
that the amplitude of longitudinal oscillations decreases. 
It would not, however, imply a diminishing SHE since it 
is related to the longitudinal momentum and not posi- 
tion. 

A closer look at the two figures reveals that in the 
transverse motions there is an additional oscillating mo- 



tion different from the spin current being in phase with 
the BOs. This is believed to be Zitterbewegung clearly 
seen in the freely evolving wavepackets of Fig. [2J Con- 
trary to that figure, here these oscillations do not collapse 
suggesting that it might be preferable to detect the Zit- 
terbewegung effect within a tilted lattice, compared to 
the case of freely propagating atoms. The survival of 
Zitterbewegung should as well derive from the alternat- 
ing motion due to the optical lattice. 

Figures [5] and [7J have been calculated for an approxi- 
mated initial ground state of the untilted lattice. As al- 
ready pointed out, experimentally the condensate might 
be prepared in the presence of an external trap, and then 
the optical lattice as well as the coupling of the internal 
states are turned on in such a way that the condensate 
is not in a true ground state of the coupled system. To 
analyze the effect of such excitations, we have redone the 
simulations of Figs. [5] and [7J for an initial Gaussian mo- 
mentum state of a form depicted in Eq. (|10[) with pa = 
and A p = 0.1. The results are very similar to the ones 
of Figs. [6] and [JJ and we therefore do not present them 
here. We have also investigated the effect of a shallow 
harmonic potential in the transverse direction during the 
evolution in the tilted lattice. Again, a transverse spin 
current is obtained, however this time the amplitude of 
transverse oscillations is decreased due to the confining 
potential. This reduction of spin separation would indeed 
be present in any experiment utilizing a confining trap- 
ping potential. The idea of employing an optical lattice 
potential as a trapping potential circumvent such short- 
ing. In order to achieve trapping in this case one must 
stay within the BO regime, i.e. adiabatic. As Figs. [6] 
and [7J demonstrates, even though the atoms oscillate a 
relatively large spin current can build up during time. 

In relation to the measurability of the spin separation, 
we may note that even though \(y)i — (y}2\ is consider- 
ably larger than the lattice spacing, it does not imply 
that the two wave functions ipi(x, y, t) and ip2(x, y, t) do 
not spatially overlap. In fact, the overlap is relatively 
large and fluorescent measurement of a single experimen- 
tal run would therefore not contain sufficient information 
to reveal the spin current. However, this is hardly never 
the case in true experimental realizations and instead an 
ensemble of measurements should be considered, where 
the ensemble average gives \{y)i — (y)z[- Since this value 
is presumably larger than the lattice spacing, it should 
indeed be measurable [6J| . In the cold atom community, 
time-of-flight measurements are often employed, in which 
it is the atomic momentum density that is being mea- 
sured. If the separation of the momentum wave functions 
is non-zero, the position space wave functions ipi(x,y,t) 
and tp2{x,y,t) will set apart further during the time-of- 
flight. If this separation is faster than the intrinsic wave 
function broadening, a finite time-of-flight measurement 
may enhance the detection efficiency. 
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FIG. 7: (Color online) Same as in Fig. [(J but with an increased 
atom-atom interaction; g = 10. 



IV. CONCLUSION 

In this paper, we have considered an ultracold atomic 
system possessing three well-studied properties of con- 
densed matter theory, Bloch oscillations, intrinsic spin 
Hall currents, and Zitterbewegung.Our motivation for 
this study comes from the fact that Bloch oscillations 
have been verified in the same type of setups [H, l3a |. 
light induced gauge potentials have been demonstrated 
[HI EH HII i and Zitterbewegung have been detected in 
ion-trap experiments [23[. 

The spin-Hall effect (SHE) and the Bloch oscillations 
(BO) in our system are driven by two perpendicular 
forces: the extrinsic force appearing due to tilting of 
the lattice potential and the intrinsic force originating 
from the effective gauge potential caused by the spin- 
orbit (SO) coupling. We demonstrated that the interplay 
between these two forces gives rise to a correlated mo- 
tion in the xy-plane, and despite the fact that the BOs 
are expected to prevent a build up of large transverse 
spin currents, we argued that the SHE should indeed be 
detectable. We utilized system parameters that should 
be achievable in current experiments, e.g., the amplitude 
of the lattice potential is five recoil energies compared to 



the experimentally utilized ones of Ref. [25| ranging from 
one to 14 recoil energies, and the SO coupling two recoil 
energies compared to 8 recoil energies in Ref. 19]. For a 
potential amplitude as large as 14E r , the dispersions are 
almost flat and longitudinal motion is almost frozen. The 
momentum along the longitudinal direction do, however, 
oscillate within the Brillouin size and, interestingly and 
somewhat surprisingly, one would therefore see an oscil- 
lating transverse spin Hall current despite lack of longi- 
tudinal currents. We furthermore showed that the effects 
are insensitive to experimentally relevant atom-atom in- 
teraction strengths. Entering a strongly interacting gas 
seemed to affect our results and we argue that this is 
likely to derive from a non-linearity induced breakdown 
of adiabaticity around the Dirac points. An additional 
interesting aspect being discussed is the analogue of rel- 
ativistic Zitterbewegung. Our numerical results indicate 
that it may be preferable to consider an optical lattice 
configuration when studying this effect since it persists 
over longer time periods. The optical lattice convey as 
well the possibility to avoid confining trapping potentials 
preventing large spin currents. 

We should also point out that the system presented in 
this paper allows for many extensions and different phe- 
nomena to be studied. For example, vortex generation 
in lattices [HI, |65j , Bose-Einstein condensation into non- 
zero momentum states [Hf, topologically ordered states 
of matter 67], or more strongly correlated systems, by 
means of second quantization of the mean-field Hamilto- 
nian, may find realization in atoms possessing an internal 
tripod structure and trapped in optical lattices. The flex- 
ibility permits as well for various extensions, for example 
time-dependent lattices or SO couplings. A recent pa- 
per considered Zitterbewegung of ultracold atoms for a 
driven SO coupling 68]. The same idea can be utilized 
in order to increase the transverse spin current in our 
model. 
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